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Abstract 

We study nonintersecting Brownian motions with two prescribed starting and end- 
ing positions, in the neighborhood of a tacnode in the time-space plane. Several expres- 
sions have been obtained in the literature for the critical correlation kernel K tac (x,y) 
that describes the microscopic behavior of the Brownian motions near the tacnode. 
One approach, due to Kuijlaars, Zhang and the author, expresses the kernel (in the 
single time case) in terms of a 4 x 4 matrix valued Riemann-Hilbert problem. Another 
approach, due to Adlcr, Ferrari, Johansson, van Moerbeke and Veto in a series of 
papers, expresses the kernel in terms of resolvents and Fredholm determinants of the 
Airy integral operator acting on a semi- infinite interval [a, oo), involving some objects 
introduced by Tracy and Widom. In this paper we prove the equivalence of both 
approaches. We also obtain a remarkable rank-2 property for the derivative of the 
tacnode kernel. 

Keywords: resolvent operator, Fredholm determinant, Tracy-Widom distribu- 
tion, Riemann-Hilbert problem, Lax pair, Painleve II equation, Brownian motion, 
determinantal point process. 

1 Introduction 

Recently several papers appeared that study n non-intersecting Brownian motion paths 
with prescribed starting positions at time t = and ending positions at time t = 1, 
see [H El SI H31 H51 UHl EH [22I EH] among many others. If n -»■ oo then the paths 
fill up a well-defined region in the time-space plane. By fine-tuning the parameters, we 
may create a situation with two groups of Brownian motions, located inside two touching 
ellipses in the time-space plane: see the third picture of Figure [TJ We are interested in the 
microscopic behavior of the paths near the touching point of the two ellipses, i.e., near the 
tacnode. 

It is well-known that the positions of the Brownian motions at any fixed time t £ (0, 1) 
form a determinantal point process. The process has a well-defined limit for n — > oo 
in a microscopic neighborhood of the tacnode. The limiting process is encoded by a 
two- variable correlation kernel K tlic (x,y) which we call the tacnode kernel. It depends 
parametrically on the scaling that we use near the tacnode. There exists also a multi-time 
extended version of the tacnode kernel [21 [3j H5| [^T] but we will not consider it in this 
paper. 
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Figure 1: n = 20 non- intersecting Brownian motions at temperature T = 1 with two 
prescribed starting and two ending positions in the case of (a) large, (b) small, and (c) 
critical separation between the endpoints. The horizontal axis stands for the time, t G 
[0,1], and the vertical axis shows the positions of the Brownian motions at time t. For 
n — > oo the Brownian motions fill a prescribed region in the time-space plane which is 
bounded by the boldface lines in the figure. In the case (c), the limiting support consists 
of two touching ellipses which touch each other at a critical point which is a tacnode. 

The tacnode kernel K taiC (x, y) can be expressed using resolvents and Fredholm determi- 
nants of the Airy integral operator acting on a semi- infinite interval [a, oo). This approach 
was followed in the symmetric case by Adler, Ferrari, van Moerbeke [2] and Johansson [21] 
and in the non-symmetric case by Ferrari and Veto [16]. Here the 'symmetric case' means 
that the two touching groups of Brownian motions have the same size, as in Figure [TJ and 
the non-symmetric case means that one group is bigger than the other. Similar methods 
were used to study a double Aztec diamond [3] . 

An alternative expression for the tacnode kernel can be obtained from the Riemann- 
Hilbert method. This approach was followed by Kuijlaars, Zhang and the author in |13j . 
In that paper we express the tacnode kernel in terms of a 4 x 4 matrix-valued Riemann- 
Hilbert problem (RH problem) M(z), which yields a new Lax pair representation for the 
Hastings-McLeod solution q(x) to the Painleve II equation. Recall that the Painleve II 
equation is the second-order, ordinary differential equation 

q"(x) = xq(x) + 2q 3 (x), (1.1) 

where the prime denotes the derivative with respect to x. The Hastings-McLeod solution 
|2"U] is the special solution q(x) of (jl.l|) that is real for real x and satisfies 

q{x) ~ Ai(x), x — > +00, (1.2) 

with Ai the Airy function. We note that the usual Riemann-Hilbert matrix *k(z) associated 
to the Painleve II equation, due to Flaschka and Newell |17] . has size 2x2 rather than 
4x4. 



2 



The RH matrix M{z) from the previous paragraph has been the topic of some re- 
cent developments. It was used to study a new critical phenomenon in the two- matrix 
model [13], and to establish a reduction from the tacnode kernel to the Pearcey kernel |19j . 
It was also extended to a hard-edge version of the tacnode [H] . 

Summarizing, there exist several, apparently different, formulas for the tacnode kernel 
K tac (x,y). It is natural to ask about the equivalence of these formulas. 

There is an interesting analogy with a model of non- intersecting Brownian excursions, 
also known as watermelons (with a wall) . The model consists of n Brownian motion paths 
on the positive half-line with a reflecting or absorbing wall at the origin. The paths are 
forced to start and end at the origin and the interest lies in the maximum position x max 
reached by the topmost path during the time interval t £ [0, 1]. 

Recently, several results were obtained about the joint distribution of the maximum 
position x max and the maximizing time £ max for such watermelons. One approach, due to 
Moreno-Quastel-Remenik |24j . involves resolvents and Fredholm determinants of the Airy 
operator acting on an interval [cr, oo). Another approach, due to Schehr [25] , involves the 
2x2 Riemann-Hilbert matrix ^(z) associated to the Hastings-McLeod solution to the 
Painleve II equation [171 [20]. The equivalence of both approaches has been established 
in a recent work of Baik-Liechty-Schehr |7j. Along the way they obtain Airy resolvent 
formulas for the entries of the RH matrix *&(z), see also [6, Sec. 1.1.3] for a similar result. 

Inspired by the work of Baik-Liechty-Schehr [7J, in this paper we obtain Airy resolvent 
formulas for the 4x4 RH matrix M(z). This is the content of Theorem 12.71 below. In 
Theorem 12.61 we will use these formulas to prove the equivalence of the tacnode kernels in 
[21 [161 [21] and [13] respectively. We also obtain a remarkable rank-2 property for the 
derivative of the tacnode kernel, see Theorems 12.31 and 12.51 

The Airy resolvent formulas will describe the entries in the first and second column of 
the RH matrix M(z), where z lies in a sector around the positive imaginary axis. We were 
unable to obtain similar formulas for the third and fourth column of M{z). The latter 
columns are relevant because they appear in a critical correlation kernel for the 2-matrix 
model [14] , It is an open problem to obtain Airy type formulas in that case. 

Another open problem is to extend our results to the hard-edge situation discussed in 
[TT] . and to extend the RH approach to the setting of the multi-time extended tacnode 
kernel. See [8] for some recent developments that could be useful in the latter direction. 

2 Statement of results 

2.1 The Riemann-Hilbert problem for M(z) 

In this section we recall the RH problem for M(z) from [131 H4j. Fix two numbers v?i,^2 
such that 

< <pi < v? 2 < vr/3. (2.1) 
Define the half-lines r&, k = 0, . . . , 9, by 

T = R+, r x = e ivi R+, T 2 = e^ 2 M + , T 3 = e i ^ ip2) R + , T 4 = e i(7r - Vl) lR + , (2.2) 
and 

T 5+k = -T k , fc = 0,...,4. (2.3) 

All rays T^, k = 0, ... ,9, are oriented towards infinity, as shown in Figure [21 We denote 
by flk the region in C that lies between the rays and r^+i, for k = 0, . . . , 9, where we 
identify T w := r . 
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Figure 2: The figure shows the jump contours in the complex z-plane and the corre- 
sponding jump matrix on Tj, k = 0, . . . , 9, in the RH problem for M = M(z). We 
denote by the region between the rays and F^+i- 



We consider the following RH problem. 

RH problem 2.1. We look for a matrix valued function M : C \ ^Ufc=o^fc) ~~ ^ C 4x4 
(which also depends on the parameters ri,r% > and si, S2) r G satisfying 

(1) M(z) is analytic (entrywise) for z G C \ ^Ufc=o^ fe ) ■ 
fi2) For z G Tfc, the limiting values 

M + (z) = lim M(z), M-(z) = lim M(z) 

z^>z,z on +-side ofY^ z^>z,z on — side ofT^ 

exist, where the +-side and —-side ofT^ are the sides which lie on the left and right 
ofTk, respectively, when traversing according to its orientation. These limiting 
values satisfy the jump relation 

M+(z) = M-(z)J k (z), fc = 0,...,9, (2.4) 

where the jump matrix Jk(z) for each ray r& is shown in Figure® 

(3) As z — > oo we have 

M(z) = (/ + ^ + ^ + O (1) ) diag((-,)-V 4 , *-V* (-^i/*, ,V 4 ) 

X^diag^V^V^V 4 ^), (2.5) 
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where the coefficient matrices M\,M2, ■ ■ ■ are independent of z, and with 
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(2.6) 



and 



e 1 {z) 

02(Z) 

e^z) 



.| ri (_ z )3/2_ 2si( _ z) l/2 +r 2 TZi 

-|r 2 2i 3 / 2 — 2s 2 z 1/ ' 2 —r%rz, 

In 

§r 2 z 3 / 2 



(_ z )3/2 +2si( _ z) l/2 



-r 2 T£. 



(2.7) 



Here we use the principal branches of the fractional powers. 

(4) M{z) is bounded as z — > 0. 

We will sometimes write M(z) = M(z; r±, r 2 , s\, s 2 , r) to indicate the dependence on 
the parameters. We could assume r 2 = 1 without loss of generality by a simple rescaling 
of z. 

The RH problem 12. II was introduced in [13] with r = in (|2.7p . The parameter r was 
introduced in |14j in the symmetric setting where r\ = ri = 1 and si = s 2 . The general 
non-symmetric case with the extra parameter r has not been considered before in the 
literature. The following result can be proved as in [13]; see also the discussion following 
Lemma 16.71 below. 

Proposition 2.2. (Solvability.) For any ri,r 2 > and si,s 2 ,r G R f/iere is a unique 
solution M(z) = M{z; r±, r 2 , si, s 2 ,t) to £/ie RH problem\2J] 



2.2 The tacnode kernel and its derivative 

Let ri,r 2 > and si,s 2 ,r Gibe fixed parameters. Let M(z) be the restriction of M{z) 
to the sector z G f2 2 around the positive imaginary axis. We extend M(z) to the whole 
complex z-plane by analytic continuation. This analytic continuation is well-defined, since 
the product J3J4 • • • J9J0J1J2 of the jump matrices in the RH problem 12.11 is the identity 
matrix. 

The tacnode kernel is defined in terms of the RH matrix M(z) by [13 Def. 2.60 

K tac (u,v) = — -4 r(0 1 l)M-\v)M(u)(l 1 0) T , (2.8) 

where the superscript T denotes the transpose. It will be convenient to denote by 

p(z)=M(z)(l 1 0) T GC 4xl (2.9) 

the sum of the first and second column of M{z). Observe that ()2.8fl — f)2.9|) both depend 
on the parameters r%, r 2 , s%, s 2 , r. 

Note that the kernel f|2.8[) has an 'integrable' form, due to the factor u — v in the 
denominator. Interestingly, this factor cancels when taking the derivative with respect to 
s± or s 2 . This is the content of the next theorem. 

*We note that [Ml Def. 2.6] has a typo: it has 1 M~ 1 («)M(d)' instead of 'M _1 (?;)M(m)'. 
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To state the theorem we parameterize 



si =: 0i s, s 2 =: 2 s, (2-10) 

where 01,02 are fixed and s is variable. In the symmetric case where s± = s 2 ='■ s, we 
could simply take 01 = 02 = 1. We also consider r\,r2 > to be fixed. Then we write 
Ktac(u, v; s, r), p(z; s, r) etc., to denote the dependence on the two parameters s and r. 

Theorem 2.3. (Derivative of the tacnode kernel.) With the parametrization (|2.10p , the 
kernel (|2,8p satisfies 

d 1 

— K tac (u,v;s,T) = (<riPi(u;s,T)pi(v; s, -r) + o- 2 p 2 (u; s, r)p 2 (v ; s, -r)) , (2.11) 

OS IT 

where pj, j = 1, ... ,4, denotes the jth entry of the vector p in (|2,9p . Consequently, if 
01 , 02 > then 

1 f°° 

K tac (u,v;s,T) = - / ((ripi(u;7,T)pi(v;7,-T) + a 2 p 2 (u; s,r)p 2 (v; s, -t)) ds. (2.12) 

Theorem 12.31 is proved in Section [3l Later we will apply the theorem with the param- 
eters 01 , 02 in (|4.39p . 

Note that the right hand side of f)2.11j) is a rank-2 kernel. We will see in Section 12.41 
below that a similar property holds when starting with the Airy type expressions for the 
tacnode kernel due to Adler et al. [21 131 HU1 12T] . 

We point out that (|2.1ip ~ ()2.12p have an analogue for the kernel K-g,(u,v) which is 
associated to the 2x2 Flaschka-Newell RH matrix *&(z), see Eq. (1.22)]. The kernel 
Kiy(u,v) occurs in Hermitian random matrix theory when the limiting eigenvalue density 
vanishes quadratically at an interior point of its support O I10j . 



2.3 Airy resolvents and the tacnode kernel 

Denote by Ai(x) the standard Airy function and by 

T ^ , , Ai(x)Ai'(y) - Ai'(x)Ai(y) f°° ... .... , , /nio , 

K Ai (x,y) = — ^ ^ ^ — y ±l= Mix + z)M(y + z) dz (2.13) 

x-y Jo 

the Airy kernel. For a € K. let 

K AUa (x,y)= Ai(x + z + <r)Ai(y + z + 0) dz (2.14) 
Jo 

be the Airy kernel shifted by 0. Let KAi,o- be the integral operator with kernel K A ^ a 
acting on the function space L 2 ([0, 00)). The action of the operator KAi, CT on the function 
/ is defined by 

poo 

[K Aita f](x)= K Aij<J (x,y)f(y)dy. 
Jo 

Define the resolvent operator RAi.o- on L 2 ([0,oo)) by 

R Ai)CT := (1 - KAi,.)- 1 - 1 = K A i, CT (l - Ka^)- 1 = (1 - Ka^)- 1 ^,,, (2.15) 

where 1 stands for the identity operator on L 2 ([0, 00)). It is known that R A i )CT is again an 
integral operator on L 2 ([0, 00)) and we denote its kernel by R a (x,y): 

/•oo 

[R A i,af}(x)= R a (x,y)f(y)dy. (2.16) 
Jo 
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We will sometimes use the notation 



(1 - Ka^)- 1 ^, V) = (1 + RavOOc, 2/) := S(x -y) + R a {x, y), (2.17) 

with 5{x — y) the Dirac delta function at x = y and R a the Airy resolvent kernel (|2.16p . We 
will often use the symmetry of the kernel, R a (x, y) = R a (y, x). Finally, we will abbreviate 
R a by R if the value of a is clear from the context. 

In a series of papers [21 [31 [161 EI]) Adler, Ferrari, Johansson, van Moerbeke and Veto 
study the tacnode problem using Airy resolvent expressions. We focus in particular on the 
paper by Ferrari- Veto [16j on the non-symmetric tacnode. Hence the two touching groups 
of Brownian motions at the tacnode are allowed to have a different size. The paper |16j uses 
a parameter A > that quantifies the amount of asymmetry, with A = 1 corresponding 
to the symmetric case treated in [21 [3], [21] . These papers also use a parameter a > that 
controls the strength of interaction between the two groups of Brownian motions near the 
tacnode. In the present paper we will denote the latter parameter with a capital E. The 
parameter S has a similar effect on the tacnode kernel as the temperature parameter used 
in [HME] (suitably rescaled). In order to be consistent with [13], we will use the notation 
a to denote 

, 7 = A 1 / 2 (l + A- 1 / 2 ) 2 / 3 S. (2.18) 

(What we call a was called a in [21 [31 [161 [21].) The papers [21 El [161 121] consider a multi- 
time extended tacnode kernel with time variables ri,T2. We will restrict ourselves to the 
single time case t\ = T2 =: r. 

With the above notations, define the functions [16J3 

b T ,z( x ) = ex P {~ T V + r3 / 3 ) Ai (?/)> witn y ■■= z + Cx + Y, + t 2 , 

b T , z (x) = exp (— VAry + Ar 3 /3) Ai(A 1 / 6 y), with y := —z + Cx + \/A(S + r 2 ), ^ 2 ' 19 ^ 



where 

C = (1 + A~ 1/2 ) 1/3 . (2.20) 

Note that in the symmetric case A = 1, we have b T)Z {x) = b T - z {x). The functions (|2.19D 
also depend on A, a (recall (|2.18p ) but we do not show this in the notation. Next, we 
define the functions 

A T , z {x) = b T , z (x) - A 1/6 / M(x + y + (r)b T)Z (y) dy 

Jo f oo (2.21) 
A T ,z( x ) = b rAx) ~ A~ 1/6 / Ai(x + y + a)b T>z {y) dy. 



Again we have A T ^ z {x) = A T - Z {x) in the symmetric case A = 1, and we suppress the 
dependence on A, a from the notation. 

We are now ready to introduce the tacnode kernel vC tac (u, v) = Cta, c (u, v; a, r) of Ferrari- 
Veto [16^1 . restricted to the single-time case t\ = T2 = r. The kernel can be represented 
in several equivalent ways. We find it convenient to use the following representation. 



^The notations b^ tT+ ^(x + a) and &^i/3 T X 2/3 a _ x i/6^( x + °") in IB] correspond to our notations 
X 1 ^ 6 b- T ^(x) and A _1//6 6„ Tj j(x) respectively. 

^his kernel is called C t £ (n, ^i, ts, ^2) in [16], with = u, £2 — v and n = T2 = r. Recall that we 
use a in a different meaning. 
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Proposition 2.4. Fix A > 0. The tacnode kernel Ct ac (u,v) of Ferrari- Veto \16\j in the 
single-time case T\ = t 2 = r can be written in the form 

C tac (u,v;a,r) = X^e^^K^ a2/3(s+t2) (-A 1 / 6 U , -X^v) 

roc f'OO 

+ (l + A ' 1 /2 ) i/3 / / ^ 1 -K AiiCT )- 1 (x,y)A T , u (x)A- T , v (y)dxdy, (2.22) 
Jo Jo 

with the notations (gjj]) and (f2TT7j) -~ (f2T2T]) . 

Proposition 12.41 is proved in Section 14.11 It was obtained in the symmetric case A = 1 
by Adler, Johansson and van Moerbeke [3j Theorem 1.2(i)]. 

In the above proposition we use the notation (|2.17p . Hence the double integral can be 
rewritten as 

poo rco 

/ / (1 - K A i >a )~ 1 (x,y)A r ,u(x)A- T}V (y)dxdy 
Jo Jo 

= R a {x,y)A T)U {x)A- T , v (y) dxdy + / A T)U (x)A- T)V (x) dx. 

Jo Jo Jo 

2.4 Derivative of the Ferrari- Veto tacnode kernel 

Inspired by Theorem 1 2. 3 1 we consider the derivative of the tacnode kernel £ta,c(u, v) with 
respect to the parameter a. 

Theorem 2.5. Fix A > 0. The tacnode kernel Ct &c {u,v) in ()2.22j) satisfies 

-^C ta ,c(u,v;a,T) = -C~ 2 (A 1/3 pi(u; a, r)pi(v ; a, -r) + X~ 1/2 p 2 (u; a, r)p 2 (v; a, -r)) 

(2.23) 

and consequently 

Ctac{n,v;a,T) = C~ 2 j (V /3 pi(ii; s, r)pi{v; s, -r) + X~ 1/2 p 2 (u; s, r)p 2 (v; s, -rfj ds. 

(2.24) 

Here we denote C = (1 + A -1 / 2 ) 1 / 3 and 

roo ^ 

pi(z;a,r)= / (1 - K A i, CT ) _1 (a;, 0) A T , z (x) dx 

J °roc (2.25) 

p 2 (z;a,r) = / (1 - K A i, CT ) _1 (a;, 0) A T , z (x) dx, 
Jo 

with the notation (|2.17p . 

Theorem 12.51 is proved in Section 14.31 In the proof we will use certain functions 
(14. 8 1) that were introduced by Tracy- Widom We will also obtain some alternative 

representations for (|2.25p . see Lemma 14.31 Note that the right hand side of (|2.23p is again 
a rank-2 kernel. This is a result of independent interest. 

By combining Theorems 12.31 and 12. 5| we obtain 

Theorem 2.6. ( Connection between tacnode kernels.) Fix A > 0. The kernel Ct &c {u, v; a, r) 
in ()2.22p equals the RH kernel Kt ac (u,v;ri,r 2 , s\, s 2 ,t) in (j2.8|) with the parameters 

n = X l '\ 7-2 = 1, Si = 1a 3 / 4 (S + t 2 ), S2 = I(S + t 2 ), (2.26) 

where we recall (|2,18D . 

Theorem 12.61 is proved in Section [4.41 
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2.5 Airy resolvent formulas for the 4x4 Riemann-Hilbert matrix 

In the proof of Theorem 12.61 we will need Airy resolvent formulas for the entries in the 
first two columns of the RH matrix M{z) = M(z; n, r 2 , si, s 2 , r). The existence of such 
formulas could be anticipated by comparing Theorems 12.31 and 12.51 

The formulas below will be stated for general values of r\, r 2 , s±, S2,t. The reader who 
is only interested in the symmetric case r\ = r 2 and s\ = S2 can skip the next paragraph 
and move directly to (|2.29p . 

For general n, r 2 > and s\, S2,t £ R, we define the constants 

C = (rr 2 + r 2 - a ) 1/8 , 

D = ,/^expf^^T 3 + 2(r 2S2 -r lSl )rV (2.27) 



and the functions 



b z (x) = V2^rl /6 exp (-r|r (z + Cx)) Ai ( r 2 2 ,Z (z + Cx + 2^ 
b z (x) = v&} /6 exp (rjV (z - Cx)) Ai (r\ ,Z f-z + Cx + 2^ 



(2.28) 



The above definitions of C, a are consistent with our earlier formulas (|2.20p and (|2.18p 
under the identification (|2.26[> . Similarly, the formulas for b z (x),b z (x) reduce to the ones 
in (|2.19p up to certain multiplicative constants (independent of z,x). 

In the symmetric case where t\ = r 2 = 1, s\ = s 2 =: s, the above definitions simplify 

to 

C = 2 1 / 3 , 

D = 1, 

a = 2 5 / 3 s - 2 2 / 3 r 2 , (2.29) 

b z (x) = \/2^exp(-r^ + 2 1/3 x)) Ai^ + 2 1/3 j; + 2s) , 

b z (x) = b- z (x). 

The fact of the matter is 

Theorem 2.7. (Airy resolvent formulas for the 4x4 RH matrix.) Denote by Mj^(z) the 
(j,k) entry of the RH matrix M(z) = M(z\ r%, r 2 , s%, s 2 , r). TTten £/te entries in the top 
left 2x2 6/oc/c 0/ M (z) can 6e expressed by the formulas 



pOO ^ 

M hl (z)= / (l-K Ai , a r 1 (x,0)b z (x)dx, 
Jo 

poo poo 

M 2 ,i(z) = -D / (l-K Ai<a )- 1 (x,0)Ai(x + y + 
Jo Jo 

poo poo 

M 1)2 (z) = -D~ l / (l-K Ai , CT )- 1 (x,0)Ai(x + y + C j)6 z (y)dxdy, 
io Jo 



a)b z (y)dx dy, 



(2.30) 



M 2 , 2 (z) 



/ (l-K Ai , (T )~ 1 (x,0)6,(x)drE, 
Jo 



where we use the notations (|2.17p and ()2.27j) — (|2.28|) (or ()2.29p in iae symmetric case). 
The entries in the bottom left 2x2 of M(z) can be obtained by combining the above 

expressions with equations (|2.36|) - (|2.37|) below. 
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Theorem 12.71 is proved in Section [5j The proof makes heavy use of the Tracy- Widom 
functions defined in Section 14.21 



Corollary 2.8. The first two entries of the vector p(z) in (|2,9p are given by 

/oo ^ 
(l-K Ai , a y\x,0)A z (x)dx, 



(2.31) 



where 



(2.32) 



P2 (z)= / (l-K Aha )- i (x,0)A z (x)dx, 
Jo 

r°° ~ 
A z {x) = b z (x) -D Ai(x + y + a)b z (y) dy, 

A z (x) = b z (x) - D- 1 / Ai(x + y + a)b z (y) dy. 
Jo 

2.6 Differential equations for the columns of M(z) 

In the proof of Theorem 12.71 we will use the differential equations for the columns of the 
RH matrix M(z) (or M(z)). Interestingly, the coefficients in these differential equations 
contain the Hastings-McLeod solution q(x) to Painleve II. We also need the associated 
Hamiltonian 

u{x) := {q'(x)) 2 - xq 2 {x) - q 4 (x). (2.33) 

Proposition 2.9. (System of differential equations.) 

(a) Let the vector m(z) = m(z;ri,r2,si,S2,T) be one of the columns of M(z), or a 
fixed linear combination of them, and denote its entries by rrij(z), j = 1, . . . , 4. Then with 
the prime denoting the derivative with respect to z, we have 

r^ 2 m'[ = 2Tmi + C 2 D~ l q{a)m' 2 

+ [Cq 2 (a) -z + 2 Sl /n - r 2 r 2 ] mi - [CZT Y(<r)] m 2 , (2.34) 



r m 



C 2 Dq(a)m' 1 - 2rm! 2 



+ [Cq 2 (a) +z + 2s 2 /r 2 - r 2 2 r 2 ] m 2 - [CDq'(a)] mi, (2.35) 
r\im^ = m'x — (C~ l u{a) — s\ + r 2 T)mi — C~ 1 D~ 1 q(a)m 2 , (2.36) 
r 2 im,i = m! 2 + C^ 1 Dq(a)mi + (C^ 1 u(a) — s 2 + r 2 T)m 2 , (2.37) 

where the constants C,D,a are defined in (|2.27p and we denote by q,u the Hastings- 
McLeod function and the associated Hamiltonian. 

(b) Conversely, any vector m(z) that solves (|2.34p - (|2.37p is a fixed (independent of z) 
linear combination of the columns of M(z). 

Proposition 12.91 is proved in Section 16.41 with the help of Lax pair calculations. There 
are similar differential equations with respect to the parameters s%, s 2 or r but they will 
not be needed. 



2.7 Outline of the paper 

The remainder of this paper is organized as follows. In Section [3] we establish Theorem 12.31 
about the derivative of the RH tacnode kernel. In Section H] we prove Proposition 12.41 
and Theorems 12.51 and 12.61 about the Ferrari- Veto tacnode kernel. In Section \5\ we prove 
Theorem 12.71 about the Airy resolvent formulas for the entries of the RH matrix M(z). 
Finally, in Section O we use Lax pair calculations to prove Proposition 12.91 
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3 Proof of Theorem 



Throughout the proof we use the parametrization Sj = o~jS with er,- fixed, j = 1, 2. We also 
assume r%,r2 > to be fixed. The RH matrix M(z) = M(z;s,r) satisfies the differential 
equation 

(3.1) 



— M(z) = V(z)M(z), 

for a certain coefficient matrix V(z) = V(z;s,r). This is described in more detail in 
Section f6. 31 At this moment, we only need to know that 



V(z) 



-2iz 



/0 0\ 



o~\ 

\0 a 2 0/ 



(3.2) 



where the matrix V(0) is independent of z. 
Lemma 3.1. We have 

M-\z;s,t) = ^ 1 M T (z;s,-r)K, 
V t (z;s,t) = -KV{z;s,-t)K- x , 

where the superscript T denotes the transpose and 



K 



I 2 
-h 



(3.3) 
(3.4) 



(3.5) 



with I2 denoting the identity matrix of size 2x2. 

Proof. The first equation follows from Lemma [6.3l below. The second equation then follows 
by combining (j3.3|) and ()3.ip . □ 

We are now ready to prove Theorem 12.31 We write (|2.8p in the form 



K tajC (u,v;s,T) 



2iti(u — v) 



^110 0) M t (v;s,-t)KM(u;s,t) (l 1 0) J 



-p T (v; s, -t)K p(u; s, t) 



2iri(u — v) 

where we used (|3,3p and the definition of p(^) (|2.9|) . From (|3.ip we then obtain 



(3.6) 



d_ 

ds 



K tac (u,v;s,r) 



2iri(u — v) 
1 

2iri(u — v) 



p T (v; s, —t) [V t (v; s, -t)K + KV(u; s, r)] p(n; s, r) 



p T (v; s, -t)K [-V(v; s, t) + V(u; s, r)] p(u; s, r) 



: P T (v;s,-t)K 




S 



p(u;s,r) 



TT 

1 T, 0\ 
= --p (u;s,-r) I Q J p(«;s,T), 

with S := diag(o"i, (T2), where the second equality follows from (|3.4p . and the third equality 
uses (|3.2|) . This proves (|2.11|) . By integrating this equality, we obtain (|2.12p . due to the 
fact that the entries of p(z;s,r) go to zero for s — > +00 if o"i,(72 > 0. The latter fact 
follows from [131 Sec. 3] if r = and is established similarly for general r£l. This ends 
the proof of Theorem 12.31 □ 
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4 Proofs of Proposition 12.41 and Theorems 12. 5L 12.6 
4.1 Proof of Proposition [2741 

Denote by A a the operator on L 2 ([0, oo)) that acts on the function / G L 2 ([0, oo)) by the 
rule 

fOO 

(4.1) 
(4.2) 



/■oo 

[A ff /](x) = / Ai(x + y + a)f(y)dy. 
Jo 



Observe that 



on account of (|2.14p . 

Now the kernel Cta, c {u,v) in [16, Eq. (5)] in the single time case t\ = T2 = r has the 
forn|§ 

C tac (u,v;a,T) = e T{v ~ u) K M ^ +T 2(u,v) 

+(' I I (l-K AUa )- 1 (x,y)[A a b^ u ](x)[A a b^ v ](y)dxdy 

1 - KAi^)" 1 ^,^) [A a b TjU ] (x)&_ Tj „(y) dxdy 
+ A 1 /« e VAr( u -,)^ A . A2/3(s+r2)( _ A i/6 U) _ AV a, ) 




/o JO 

-X l ' & C 



Jo 



-A 1 / 3 C 




^l-KAi,.)" 1 ^^) 



'0 JO 

-oo roo 



A ft 



(x) 



A /) 

- t ^-a u —r,v 



(y) dx dy 



o JO 



(l-KAi,<j) H 2 ^) A CT 6 T)W (x)&_ Tj „(y)dxdy, (4.3) 



where again C = (1 + A 1 / 2 ) 1 / 3 and (|2.17p . The first term in the right hand side of (|4.3D 
equals 

poo 

e T ^K Ai ^ +T 2(u,v) = C & T , u (x)6_ Ti „(x)dx, (4.4) 

Jo 

on account of (|2.19p . The second term in the right hand side of (j4.3|) equals 

noo 
(1 - K A i, CT ) _1 (x, y) [A a b T)U ] (x) [A a b- T>v ] (y) dx dy 

noo /*oo 
(1 - KAi, CT ) -1 (x,y) b T)U (x)b^ T ^ v (y)dxdy - C \ b T , u (x)b- T)V (x) dx, (4.5) 
Jo 

where we used that 

A a (l - Ka^)- 1 A a = (1 - KAi,,)- 1 ^,, = (1 - KAi,.)- 1 - 1, 
on account of (|4.2p . Finally, the third term in the right hand side of (|4.3p can be written 



as 



X^C 




1 - K A i, CT ) (x, y) [A a b TjU ] {x)b- TyV (y) dx dy 



o Jo 



-X l ' G C 



o Jo 



(1 - K A i l(T ) 1 (x,y) b TjU (x) 



A h 



(y)dxdy, (4.6) 



§ The notations a, a, b^ +6 (x + a), B^ a+( (x + a), f^i/a^/a^i/e^ + 5), -B^i/3 TiA 2/3 CT _ A i/6 5 (a; + 
a) in 16 correspond to our notations E, a, X 1 ^ (i b^ T ^(x), [Aa-b- T ^] (x), X~ 1 ^ (i b- T ^(x), ^Aa-b- T ^~^ (x), 
respectively. 
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since the operators (1 — Kju i(T ) 1 and A CT commute. Inserting (|4.4p - (|4.6p for the first three 
terms in the right hand side of ()4.3|) . we obtain 



C ta , c (u,v;a,T) 

A l/6 e v^r( U -,)^ A . A2/3(s ^ 2)( _ A l/6 U) _ A l/6 u) 
roo poo 

+C / (1 - K A i >(T ) _1 (2;,y) b T , u {x)b- T , v {y) dxdy 
Jo Jo 



-X^C 
-X^C 



(1 - K A i, CT ) 1 {x,y) b T ^ u (x) 



o jo 

oo poo 



A /) 



(y)dx dy 



o JO 
oo /-oo 



(l-KAi,<j) 1 (x,y) A a b TiU (x) 



A h 



(y)dx dy 



o Jo 



:i-K AijCT )- 1 (x, 2/ ) 



A ft 



(x)b- T , v (y) dx dy, 



which is equivalent to the desired result (|2,22D , 



(4.7) 
□ 



4.2 Tracy- Widom functions and their properties 

In some of the remaining proofs we need the functions 

poo 

Q(x)= / (l-K AiiCT )- 1 (x, 2 /)Ai(y + a)dy 

J foc (4.8) 

p(x) = (1 - Ka^)- 1 ^, vW(v + dy, 

JO 

with the usual notation (j2. 17j) . Note that Q, P depend on a but we omit this dependence 
from the notation. The functions Q, P originate from the seminal paper of Tracy- Widom 
|26| . Below we list some of their properties. These properties can all be found in [26J, 
see also [5j Sec. 3.8] for a text book treatment. One should take into account that our 
function R(x,y) = R a (x,y) equals the one in [26], [5j Sec. 3.8] at the shifted arguments 
x + a, y + a, and similarly for Q(x) and P(x). 



+ )R( x ,y) = R(x,0)R(0,y)-Q(x)Q(y), (4.9) 



Lemma 4.1. The functions Q, P and R = R a satisfy the differential equations 

' d_ + d_ 

dx dy 

^R(x,y) = -Q(x)Q(y), (4.10) 

Q'(x) = P(x) + qR(x,0)-uQ(x), (4.11) 

P\x) = (x + a-2v)Q(x)+pR(x,0) + uP(x), (4.12) 



where 



q = Q(0) (4.13) 
p = P(0) (4.14) 

Q(x)Ai(x + a) dx (4.15) 



/*oo roc 

/ Q(x)Ai'(x + a)dx = P(x)Ai(x + a) dx. (4.16) 
Jo Jo 
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Note that q,p, u, v are all functions of a, although we do not show this in the notation. 



They satisfy the following different 



al equations with respect to a [26], Sec. 3.8] 

= p — qu 
= aq + pu — 2qv 

= V 

= -VI- 



(4.17) 
(4.18) 
(4.19) 
(4.20) 



It is known that q = q(cr) is the Hastings-McLeod solution to the Painleve II equation 
(ll.ip ~ (jl.2p . Moreover, u = u(o~) is the Hamiltonian (|2.33p . and 



2v 



2 2 

u — q . 



(4.21) 



Finally, we establish the following lemma. 
Lemma 4.2. For any function b on [0, oo) we have 



/ Q(x)b(x)dx= / / (l-K Ai)(7 )-\x,0)Ai(x + y + o-)b(y)dxdy (4.22) 
Jo Jo Jo 



and 



o Jo 



Q (x) Ai(x + y + a)b(y) dx dy = / R(x, 0)6(x) dx, 



(4.23) 



with the usual notations R = R a and (j2. 17[) . 



Proof. Recall the operator A a defined in (|4.ip ^ ()4.2p and let 5q denote the Dirac delta 
function at 0. We have 

/ Q(x)b(x)dx = / / (l-KAi^r 1 ^,?/) [A a 5 }(y)b(x)dxdy 
Jo Jo Jo 

/■oo roo 

= / / (l-K Ai , a )- 1 (x,y)5 (y) [A a b](x)dxdy, 
Jo Jo 

where we used that the operators A a and KAi l(T commute. This proves (|4.22p . Next, 

Q(x)Ai(x + y + a)b(y) dx dy 
(l-KAi,.)- 1 ^,^ [A a S ] (y) [A a b](x)dxdy 
R(x,y)d~ {y)b(x)dx dy, 

Jo Jo 
where we used that 



o Jo 



JO 

OO POO 



A CT (1 - KAi, ff ) _1 A ff = (1 - Ka^)- 1 ^^ = R Ai)CT . 



This proves ()4T23l . 



□ 
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4.3 Proof of Theorem 1231 

Lemma 4.3. The formulas (|2.25|) can be written in the equivalent ways 



pi(z;a,r) = A T , Z (P) + / R(x,0)A T , z (x) dx 

Jo 

= b T ,z(0) - A~ 1//6 / Q(x)A T , z (x)dx 
Jo 

/•oo 

p 2 (z;a,T) = Ar,z(°)+ / R(x,0)A TtZ (x)dx 

Jo 

roo ____ 
= brM-\ 1/% / Q(x)A T , z (x)dx, 



with the notations R = R a and 
Proof. Immediate from Lemma 14.21 and the definitions. 
Lemma 4.4. We have 



ox 



&r,«(ic), 



(l + A- 1 ^)^.^. (a; ) = a-V2|_^ ( a; ) + A 1 /6 Ai ( a; + CT )& (0). 



(4.24) 
(4.25) 
(4.26) 
(4.27) 

□ 

(4.28) 
(4.29) 
(4.30) 



Proof. The first two formulas are obvious from the definitions (|2.18|) - (|2.2U|) . For the last 
formula, we calculate 



1/2- 


d 


-\- 


-i/2 A" 




do 


ax 


1/2- 


d 


- A" 


-i/2 






dx 



3 -A-*l 



At.z (x) 

M*) - A 1 / 6 Ai(x + y + a)b T , z (y) dy^j 
Ai(x + y + cr)b T Jy) dy 



Ai'(x + y + CT)6 TjZ (y) + Ai(x + y + o>' (y) ) dy 



= -A 1 / 6 

= -A 1 / 6 

= A 1 / 6 Ai(x + a)6 T , z (0), 

where the second and third equalities use (|4.28p and (|4.29p respectively, and the last 
equality follows from integration by parts. This proves the lemma. □ 



Now we prove Theorem 12.51 With the help of (|4.30p , the derivative of (|2.22|) with 
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respect to a becomes 







( 1 + X-V 2 ) 2 /3—£ tac ( u , V ;a,T) 
= (1 + \-^fl*\^e^-v)^K M> x2/ s { x +t2) (-\ 1/6 u, -\V*v) 

noo 
(I-Km^-^x^) (A' T Jx)A- T> v(y) + A T , u (x)AL TtV (y)) dxdy 

+A 1/6 6 TiU (0) / Q(y)A-r, v (y)dy 
Jo 

+A 1 / 6 6_ T ,„(0) / Q(x)A T>u (x)dx 
Jo 

i — R(x, y) Ar, u (x)A- T , v (y) dx dy, (4.31) 

where we used the definition of Q in (|4.8p . The first term in the right hand side of (|4,31D 
can be written as 

(1 + A-^)2/3 A i/6 e ^( u -,)^^ A . A2/3(s+r2)( _ A i/6 U) _ x i/e v) 

= -A 1/3 e v/Xr(n - ,;) Ai(-A 1/6 u + A 2/3 (£ + r 2 ))Ai(-A 1/6 w + A 2 / 3 (£ + r 2 )) 

= -X 1/S b T J0)b- T J0), (4.32) 



where the first equality follows from (|2.18p and (|2.14|) and the last equality uses (|2.19p . 
The second term in the right hand side of (|4.3ip can be written as 



A" 1 / 2 / / (I-Km^-^x^) (A' TiU (x)A- T , v (y)+A T , u (x)A'_ T!V (y)) dxdy 
Jo Jo 

Tx + Wy) R ^ y) 



. A -i/2 




J 



A T ,u(x)A- T ,v(y) dx dy + A T , u (0)A- T ,v(0) 



+A T , u (0)J^ R(0,y)A-T,v(y) d y + A-T,v(°) J R(x,0)A T)U (x)dx) , (4.33) 
where we used ()2.17p and integration by parts. Finally, we observe that 



d ( d d 

^ do ^ \dx dy 



R(x,y) = -\-V 2 R(x,0)R(0,y) - Q(x)Q(y), 

(4.34) 



on account of (1P1)- (H3UD . Inserting (|Q2l) - (Q4"l) in the right hand side of (I4T3TD . we 
obtain 

{1 + X -V 2 f 3 -?-£ tac ( u ,v,o,T) 



- -A 1/3 (b TyU (0) - A" 1 / 6 J Q(x)A T ,u(x) dx) f 6- r , v (0) - A" 1 / 6 J Q(x)A- T , v (x) dx 

-A~ 1/2 (a t ,u(0) + J R(x, 0)Ar, u (x) dx^j (a-t, v (0) + 22(0, x)A- T)V {x) dx\ 
= -A 1/3 pi(u;o-,r)pi(t;;cr, -r) - X~ 1/2 p2(u;a,T)p 2 (v;a, -r), (4.35) 
on account of (|Q5jl - ([OBp . This proves (jT25jl . 
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Finally we prove (|2.24p . From [5, Sec. 3.8] we have the estimate 

\R*{x,y)\ <C e- x -y- 2 ° (4.36) 

for all x,y,o~ > 0, with Co > a certain constant. We also have the asymptotics for the 
Airy function 

Ai(x) ~ exp (~^ 3/2 ) /(2^x 1/4 ), x -> +oo. (4.37) 

Consequently the kernel £ tac (-u, v; a, r) in (|2.22 j) . ()2.17|) goes to zero (at a very fast rate) 
for a — > +oo. Integration of (|2.23p then yields (|2.24p . This proves Theorem 12.51 □ 



4.4 Proof of Theorem [Ml 

Let the parameters ri,r2,Si,S2 be given by (|2.26p . As already observed, the expressions 
for C, a in (|2.27p and (|2.20p . f|2. 18[) are equal under this identification. Similarly, the 
expressions for p\ , p 2 in (|2.3ip and p\ , p 2 in (|2.25p are related by 

Pj (z) = V^r 1 / 6 exp (V 4 r(£ + Hr 2 )J pj(z), j = 1, 2. (4.38) 

Note that the exponential factor in (|4.38[) cancels out in (|2.11|) . 

Now observe that the formulas for s%, s 2 in (|2.26|) are of the form Sj = OjS with 

s := (S + r 2 )/2, (7i := A 3 / 4 , a 2 := 1. (4.39) 

By (|2.18p we then have 

— Aac(«, «; ff, t) = 2\/A(l + A" 1 / 2 ) 2 / 3 — £ tac (n, «; a, r), (4.40) 
as oa 

where we consider r to be fixed. Theorem 12.61 follows from Theorems 12.31 and 12.51 and the 
above observations. □ 



5 Proof of Theorem 12.7 
5.1 Preparations 

Lemma 5.1. The functions b z (x), b z (x) in (|2.28p satisfy the differential equations 

b' z (x) = C-^b z (x), b' z (x) = -C-^b z (x), (5.1) 

and ^ 

r 2 2 ^b z (x) + ^b z (x) = (z + Cx + 2^ - r 2 r 2 ) 6,(x), (5.2) 

r ^l& hz{x) ~ 2r ^ (x) = {~ z + Cx + 2 n ~ r ' r2 ) ~ bz{x) - (5 ' 3) 

Proof. Equation (|5.1|) is obvious. The other equations follow from a straightforward cal- 
culation using the Airy differential equation Ai (x) = xAi(aj). □ 
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Lemma 5.2. The function A z {x) in (|2,32p satisfies the differential equations 





3z A z {x) = C- 1 (A' z (x) - DAi(x + a)b z (0)) (5-4) 



and 



r 2 2 ^A z (x) + 2t-^A s (x) = (z + Cx + 2^- r 2 r 2 ) A z (x) 

+ CD (Ai(x + a)b' z (0) - Ai'(x + cj)^(O)) . (5.5) 

Proof. Equation (|5.4p follows from the definition of A z , (|5.ip and integration by parts. 
Now we check the formula (|5.5p . From (|2.32p we have 

-2 d 2 

d 2 , , , , „ Z" 00 . . , , d 2 



a f°° a ~ 

= r 2Q-2 b z( x )- r 2 D j o ^-{x + y + a)^b z {y)dy 

Q2 roo q2 ^ roo ^ 

= r 2 2 in h ^ x ) + r i 2D Ai(x + y + a)—^b z (y)dy-CD Ai(x + y + a)b" z {y) dy 
dz l J dz 2 J 

where in the second equality we used rj~ 2 + r 2 2 = C 3 and (|5.ip . Hence 



d 2 . , , „ d . , , / _ 9 a 2 , , s d 



r 2 1T2 A ^ X ) + 2r— = r 2 2 7r^& 2 (a0 + 2r— b z (x] 
az 1 az V az" 1 



+ £>y Ai(x + y + a) l r - 2 —b z (y)-2r— b z (y)j dy-CDj M{x + y + o)b" z {y) dy. 

(5.6) 



In the first two terms in the right hand side of (|5.6p we use the differential equations (|5.2p ~ 
(|5.3|) . and in the third term we integrate by parts twice and subsequently use the Airy 
differential equation Ai"(x) = xAi(x). The lemma then follows from a straightforward 
calculation, taking into account that 

-z + Cy + 2— - r 2 r 2 ) - C (x + y + a) = -z - Cx - 2— + r 2 2 T 2 (5.7) 
n J r 2 

thanks to (j23T)) . □ 

Lemma 5.3. The formulas (|2.3ip can be written in the equivalent ways 

rOO 

Pi(z)=A z (0)+ R(x,0)A z (x)dx (5.8) 
Jo 

= b z (0)-D- 1 Q(x)A z (x)dx, (5.9) 



roo 

P2(z)=A z (0)+ R(x,0)A z (x)dx (5.10) 

J 

roo ____ 

= b z (0)-D Q(x)A z (x)dx. (5.11) 
J o 

Proof. Immediate from Lemma 14.21 and the definitions. □ 
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5.2 Differential equation for pj 

In this section we check that the expressions for p\(z), p%(z) in (j5.9j) (j5. 10p satisfy the 
differential equation (|2.34p (with rrij :=pj). From (|5,9p we have 



P'l(z) = §~ z ["6,(0) 



D~ 



Q(x)— [A z (x)] dx. 



(5.12) 



We calculate the second term 







Q{x)— [A z {x)\ dx 



C- 1 I J Q(x)A' z (x) dx - Dub z (0) 
-C- 1 (J Q'(x)A z (x) dx + qA-M + Dub z (0) 

-C- 1 (J [P(x) + qR(x, 0) - uQ{x)] A z {x) dx + qA-M + Dub z (0) 
-C' 1 (^Dupi(z) + qp 2 (z) + J P{x)A z {x) dx 



(5.13) 



where the first equality uses (|5.4p and (J4.15P , the second one uses integration by parts and 
(|4.13p . the third one uses (|4.1ip . and the fourth equality uses (|5,9p - (|5,10D , 
From (f5TT2l) - (l5TT3]l and r\ 2 + r 2 2 = C 3 , we get 



rrVi(^) = t V^ (0)] + r ^ D ' x i Q{x) d~ z [Mx)] dx 

+ C 2 D- 1 \ Dupi{z) + qp 2 (z) + J P{x)A z {x) dx 



(5.14) 



By differentiating (|5.14p with respect to z, we get 

dz 2 



r^p'l(z) - C 2 D-\p' 2 {z) - 2rp' 1 {z) = rf 2 ^ \b z (0) 



2 — 

dz 



&*(o) 



+ D- 1 Q(x) (V [M*)\ + ^Yz [Az{x)] ) dx 

C 2 D- X \Dup' x {z) + P(x)^- [A z (x)} dx^j , 
where the last term in the left hand side was expanded using (|5.12p . Equivalently, 



r^p'i(z) - C 2 D-\p' 2 {z) - 2rp' 1 (z) = r^ 2 & ~ 



dz 2 . 



6,(0) 



2r 



d_ 

dz 



6 2 (0) 



+ CD' 1 C 



Y-l 



" ( r 2 2 ^2 IM*)] + [A z (x)] ) «Lr 



[■CO Q \ 

+CDup' 1 {z) + C J P(x)^ [A z (x)} dx) . (5.15) 



We will calculate each of the terms in the right hand side of (|5.15p . We start with 



(' I P(x)^- z [A z (x)]dx 



P(x)A' z (x)dx - Dvb z {0) 

P'(x)A z (x)dx +pA z (0) + Dvb z (0) 
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where the first equality follows from (|5.4[) and (|4,16p . and the second equality uses inte- 
gration by parts and (|4.14p . Consequently, 

roc Q roc 

C / P{x)— [A z (x)] dx = - (x + a)Q{x)A z {x) dx 
Jo oz Jq 

P{x)A z {x) dx-pp 2 {z) - 2Dvp 1 (z) + Dvb z (0) (5.16) 



by virtue of g32D and 

Next we calculate the third term in the right hand side of (|5.15p . 

C- 1 Q(x) [r 2 2 ^ \A z (x)\ + 2r^- [A z (x)]\ dx 

z + Cx + 2—- rfr 2 ) Q(x)A z (x) dx + D(ub' z (0) - vb z (0)), (5.17) 



C 



on account of (|5.5|) and (|4.15|) - ([4.16|) . By adding (|5.16|) - (|5.17|) and canceling terms we 
get 

POO £3 /*00 / £l2 CJ \ 

C J P(x)— [A z (x)] dx + C- 1 J Q(x) [r 2 2 ^ [A z (x)} + 2r— [A z (x)]j dx 

Q(x)A z (x) dx 



P(x)A z (x) dx + CT 1 ( z - 2— + r\r 2 
V n 

+ Dub'M - 2Dv Pl (z) -pp 2 (z), (5.18) 

where the factor between brackets in front of the second integral in the right hand side 
was obtained via (|5.7|) . Finally, we calculate the fourth term in the right hand side of 
(l5TT5j) . 

d ~ f°° 
CDup[(z) = CDu—[b z {0)} + Du 2 pi(z) + uqp 2 (z) +u P(x)A z (x) dx 

~ r°° 
= -Dub' z (0) + Du 2 Pl (z) + uqp 2 (z) +u P{x)A z {x) dx (5.19) 



where the first equality follows from ()5.12[) ()5.13|) and the second one from (|5.ip . Inserting 
(|5.18|) - (|5.19p in the right hand side of (|5.15p and canceling terms we get 



iTVW - C 2 D' l qp' 2 (z) - 2rp' l {z) 



d 2 [~ _l „ d 



6,(0) 



^ T dz 



1 dz 2 

+C[u 2 - 2v] Pl (z) + CD- l [uq-p]p 2 {z) 



1,(0) ) - r D- 1 (z-2^-+r 2 T 2 



Q{x)A z {x) dx 



-z + 2—- r\T 2 + C(u 2 - 2v) 



. z + 2 £l - r 2 r 2 + Cg 2 
ri 



pi{z) + CD 1 [uq-p]p 2 (z) 
p 1 (z)-[CD- 1 q'] p 2 (z) 



where the second equality follows from (|5.9p and (|5.3p . and the last equality uses (|4.17p 
and (I4.2ip . We have established the desired differential equation ()2.34p . 
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5.3 Other differential equations 

Denote by Nj t k(z) the right hand side of the formula for Mjk(z) in (|2,30p . Let the vectors 
p(z),m(z) have entries 

Pj {z) = N jA {z) + N ji2 (z), rrijiz) = N^z) - N jl2 (z), 

for j = 1,...,4. For p(z) this is compatible with (|2.3ip . We have already shown in 
Section [5.21 that p(z) satisfies the differential equation (|2.34p (with rrij := pj). We claim 
that the same statement holds for m(z). This can be shown by going through the proofs 
in Sections I5.1H5.2I again and replacing the appropriate plus signs by minus signs and 
vice versa. We leave this to the reader. Summarizing, both p(z) and m(z) satisfy the 
differential equation (|2.34p . By symmetry they also satisfy the differential equation (|2.35p . 
Finally, (|2.36p - (|2,37p are valid by construction. Proposition 12. 9f b) implies that p(z) and 
m(z) are fixed linear combinations of the columns of M{z). By linearity, the same holds 
for the Nj t k(z). 

5.4 Asymptotics 

In view of the previous section, Theorem 12.71 will be proved if we can show that the 
expressions for Mj t \(z) and Mj ;2 (z) in (|2.30p and (|2.36p - (|2.37p satisfy the asymptotics for 
z — > oo in the RH problem 12.11 It will be enough to prove the asymptotics for the second 
column Mji{z). Moreover, it is sufficient to let z go to infinity along the positive real line. 
Indeed, first observe that the second columns of M(z) and M(z) are equal if Rez > 0. 
Furthermore, the second column of M(z) is recessive with respect to the other columns as 
z — > +oo, due to (|2.5p . So if we can prove that the expressions for Mj^(z) in ()2.30p share 
this same recessive asymptotic behavior for z — > +oo, then we are done. 

The proof of the above asymptotics is now an easy consequence of ()4.36p - (|4.37p . This 
ends the proof of Theorem 12.71 



6 Proof of Proposition 12.9 

6.1 Painleve formulas for the residue matrix Mi 

Let the parameters r±,r 2 > and s±, s 2 , r € R be fixed. In the proof of Proposition 12 . 91 we 
will need Painleve formulas for the entries of the 'residue' matrix M\ = Mi(ri, si, s 2 , t) 
in (|2.5|) . Write this matrix in entrywise form as 



Mi 



fa b 

-b —a 

ie if 

\if ie 



ic id\ 

id ic 

-a P 

-P 5/ 



(6.1) 



for certain numbers a,a,b,b,c,c, . . . that depend on r±, r 2 , si, s 2 , r. We will sometimes 
write a(r\, r 2 , s\, s 2 ,t), b(ri,r 2 ,si,s 2 ,r), etc., to denote the dependence on the parame- 
ters. 

In the symmetric case n = r 2 , s± = s 2 we are allowed to drop all the tildes from 
(|6.1|) . while in the case r = we can replace all the Greek letters in (|6.ip by their Roman 
counterparts and put d = d and / = /• These are special instances of the next lemma. 
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Lemma 6.1. (Symmetry relations.) Let ri,r 2 > and si,S2,t 6 R be fixed. The 16 
entries x = x{r\, r 2 , s\, s 2 , t) of the matrix (|6.ip are all real-valued and they satisfy the 
symmetry relations 

x(n,r 2 , sx,s 2 , t) = x(r 2 ,r 1 ,s 2 , si, r), (6.2) 
for any x = a,b, c, d, e, f, a, j3, and 

x{n,r 2 , si, s 2 , t) = x(n,r 2 , s 1 ,s 2 , -t), (6.3) 

/or any x = a,b, a, b, c, c, d, e, e, /, where we write x = a, /3, a, /3, c, c, d, e, e, /, respectively. 

The proof of Lemma 16. II follows from Section [6.21 below. 

Now we relate the entries of the matrix Mi to the Hastings-McLeod solution q(x) of 
Painleve II and the Hamiltonian u(x) in (|2.33p . The next theorem was proved for the 
special case r = in |13| and in the symmetric setting r\ = r 2 , s\ = S2 in 114], In the 
general case we have the extra exponential factor D in (|2.27|) . 

Theorem 6.2. Let the parameters ri,r 2 > and Si,S2,t 6 R be fixed. The entries in the 
top right 2x2 block of (|6.ip are given by 

d = {r2CD)- 1 q{a) (6.4) 

d = ( ri C)- l Dq(o-) (6.5) 

c = r^ l {s\-C- l u{a)) (6.6) 

c = ^(sl-C'M*)) (6.7) 

where q is the Hastings-McLeod solution to Painleve II (|l.ip - (jl.2p . u is £/te Hamiltonian 
(|2.33p . and with the constants C,D,a given by (|2.27[) . Moreover, some of the other entries 
in (|6.ip are given by 

b = (c + ^d-irj^D)- 1 ^^) (6.8) 
6 = ( C + Tn)d- {rlC 2 )- l Dq'{o-) (6.9) 
/3 = (c- T r 2 )d- (nr 2 C 2 )- 1 Dq'((T) (6.10) 
/3 = ( c -Tn)d-(rir 2 C 2 L>)-Y(o-) (6.11) 



and 



r 2 <9d 
r 2 + r| 9t 



J"i / = — ^ 2— + (— ric — r 2 c + r 2 T + r 2 r)o — rid 2 d + r2C 2 d — 2s 2 d (6.12) 



f2f = o ri 9 ^- + (~nc - f2C + t*?t + rlr)b - r 2 d 2 d + ric 2 d - 2s\d. (6.13) 

r^ irjOT 

Theorem 16.21 is proved in Section | 



6.2 Symmetry relations 

For further use, we collect some elementary results concerning symmetry. 

Lemma 6.3. (Symmetry relations.) For any fixed r%, r2, s\, s 2 ,t, the RH matrix M sat- 
isfies the symmetry relations 

h 0\,,, Jh 



M(z;ri,r 2 ,si,s 2 ,r) = ( Q ^ \ M(z; n, r 2 , si, s 2 , r) f Q _ J ) , (6.14) 
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where the bar denotes complex conjugation, 

M~ T (z;ri,r 2 ,si,s 2 ,T) = K~ 1 M(z;ri,r 2 , si, s 2 , -t)K, (6.15) 
where the superscript ~ T denotes the inverse transpose, and finally 

M(-z;n,r 2 ,si,s 2 ,T) = ( J Q ^ J jM(z;r 2 ,r 1 ,s 2 ,s 1 ,T)('^ _°^J , (6.16) 



where we denote 



Proof. This follows as in [1 3|. Sec. 5.1]. One easily checks that the left and right hand 
sides of ()6.14p satisfy the same RH problem. Then (|6.14p follows from the uniqueness of 
the solution to this RH problem. The same argument applies to (|6.15p and (|6.16p . □ 

Corollary 6.4. For any fixed ri,r 2 , si, s 2 ,t, the residue matrix M\ in (|2.5|) . (|6.1|) satisfies 
the symmetry relations 

M 1 (r 1 ,r 2 ,s 1 ,s 2 ,r) = ^ M x (n, r^, si, s^, r) ^ 2 , (6.18) 

M?(n,r 2 , Sl ,s 2 ,T) = -K- 1 M 1 (r 1 ,r 2 ,8 1 ,8 2 ,-T)K, (6.19) 
M 1 (r 1 ,r 2 ,s 1 ,s 2 ,r) = M x (r 2 , n, s 2 , si, r) ({ , (6.20) 

mi/t the notations J, K in (|6.17p . 

Lemma 16.11 is an immediate consequence of Corollary [6? 



6.3 Lax system 

To the RH matrix M(z) there is associated a Lax system of differential equations 

o o o 

— M = UM, —M = VM, —M = WM, (6.21) 

OZ OS OT 

for certain coefficient matrices U,V,W. These matrices were obtained in the symmetric 
case r\ = r 2 , s\ = s 2 in [Til Sec. 5.3], [H] and for r = in [131 Sec. 5.2]. We will consider 
the general nonsymmetric case. To take derivatives with respect to s\ or s 2 , we again 
parameterize Sj = ajs with a±,a 2 fixed and s variable, as in (|2.10p . 

Lemma 6.5. In the general nonsymmetric setting, with the parametrization s\ = o~\s, 
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S2 = the coefficient matrices U, V, W in (|6,2ip take the form 



U 



' —ric + rfr 

—rid 

(r\c 2 — r2dd — 2s\ + r\z)i 
-(ri/3 + r 2 b)i 



r 2 d 
r 2 c — r\r 
-(rib + r 2 /3)i 
(r2C 2 — r\dd — 2s2 — r2z)i 



r\% 


r\c + rfr 
-r 2 d 



V 



( 


0\C 


a 2 d 


—o\i 


\ 




o\d 


o- 2 c 





a 2 i 




-c 2 + ^dd+f iS -z)i 


(a\b - a 2 P)i 


— 0\C 


—a\d 


V 


(ax/3 - a 2 b)i 


a 2 ^-^dd-^s-z)i 


-a 2 d 


-a 2 cj 





r 2 i 
r\d 
-r 2 c — r\r ) 

(6.22 



(6.23) 



/ d z 

r'f+r'i 



W = (r 2 1+ rl) 
with the notations in (16.1 



~7~, — 7 > 



-fi 



di 



-di 



-P 



(6.24) 



Proof. This is a routine calculation, that follows similarly as in the above cited references 
[TT| IT3l fT4"] . Prom the asymptotics (|2.5[) we have for z — > oo that 



9 M ,r 1 

oz 



( 



1+ + ••• 



rfr 








— rST 



z(TiZ - s ly 




z(ri — si£ 

rfr 



- 1 ) o 

i{r 2 + 1 


— r2r 



-i(r 2 z + s 2 ) 

/-^ + ...^ + 0(z~ 1 ). (6.25) 



Since the RH matrix M(z) has constant jumps, the left hand side of (|6.25p is an entire 
function of z. Liouville's theorem implies that it is a polynomial in z. Collecting the 
polynomial terms in the right hand side of (|6.25)) we obtain 



U 



( r h 





ir\ 


\ 





-t\t 





ir 2 


i(r x z - si) 





rfr 





V o 


-i(r 2 z + s 2 ) 





-t\t) 



+iM 1 A-iAM 1 , A := 



/0 0\ 



ri 

\0 -r 2 0/ 



With the help of f)6. lj) and a small calculation, we then get (|6.22p . To obtain the (3,1) 
and (4,2) entries of (|6.22|) . we also need the relations 



a + a 
a + a 



-c 2 + ^dd + * 



H 

r2 



-c H aa H , 



£2 

T2 



(6.26) 
(6.27) 



which follow from the fact that the (1,3) and (2, 4) entries in the z 1 coefficient in (|6.25D 
are equal to zero. A similar argument yields ()6.23p - (|6.24p . □ 
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6.4 Proof of Proposition 12.91 

Let the vector 111(2) be a solution of J^m = Um with U in ()6,22p . By splitting this 
equation in 2 x 2 blocks we get 

P £ m #*n = f" C + r 5 ~ rV2d> ) Mil +< Mf) > (6-28) 
r 2 m 2\ z )) \-rir 2 d c - r 2 T J \m 2 (z) ) \m 4 (z)J 

r± m' 3 (z)\ . I c 2 — r^ 1 r 2 dd — 2rj~ 1 s\ + z —b — r^ 1 r 2 f3 \ /m\{z) 

r 2 1 m' 4 (z)) \ - rir -lp-l, c 2 -r 1 r 2 1 dd-2r 2 1 s 2 - zj \m 2 (z) 

+ ( C + Y J V m flV (6.29) 
V -d -c-r 2 rj \m 4 (z)J 

From (|6.28|) and (|6.4|) — (|6.7|) . we easily get (|2.36p - (j2.37p . To prove the two remaining 
differential equations, we take the derivative of (|6.28p and use (|6.29p to get 

r^ 2 m'((z)\ _ f-c + riT r^ 2 r^d\ fr^m'^zf 
r 2 2 m 2 (z)) ~ \-r 2 r 2 2 d c-t 2 t) \r 2 1 m' 2 (z) / 

(c 2 — r^ 1 r 2 dd — 2r^ 1 si + z —b — r^ l r 2 j3 \ (m\{z) 

-r\r 2 x fi-% c 2 - rir 2 1 dd-2r 2 1 s 2 - z) \jn 2 (z) 



+ 



c + r\T d 
—d —c — r 2 r 



r 1 1 m' 1 (z)\ f—c + r\T r 1 1 r 2 d\ ( m\{z 



r 2 1 m' 2 (z) J \—rir 2 l d c — r 2 T J \m 2 (z] 



From this equation and (|6,4p - (j6.1ip we obtain the desired differential equations (|2.34p ~ 
()2.35p . This proves Proposition 12.9( a) . 

To prove Part (b), let m(z) satisfy the differential equations (|2.34p ~ (|2.37p . From the 
proof of Part (a) above, we see that J|m = Um with U in ()6.22p . But then 

a 

— [M _1 m] = M _1 (C7 - U)m = 0, 
oz 

which implies Proposition 12.9( b). □ 



6.5 Proof of Theorem 16.21 

The matrices U, V, W in (|6.2ip satisfy the compatibility conditions 



dU dV , 
— = — -UV + VU 6.30 
OS oz 

dU dW 

-77- = -7T- ~ UW + WU. (6.31) 

or oz 

dzds dsdz ' 



These relations are obtained by calculating the mixed derivatives 755^ M = jf^M and 



M = a a M, respectively, in two different ways. 



dzdr drd 



Lemma 6.6. Consider the matrices U,V,W in (|6.22|) ()6.24j) . Then with the expressions 
(|6.4|) - (|6.13p in Theorem 1 6. "A the compatibility conditions (|6.3Up - (|6.31|) are satisfied. 
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Proof. This is a lengthy but direct calculation. It is best performed with the help of a 
symbolic computer program such as Maple. First we consider the compatibility condition 
with respect to s. Writing the matrix equation (|6.30p in entrywise form, with the help of 
Maple, we obtain the system of equations, with the prime denoting the derivative with 
respect to s, 

(6.32) 
(6.33) 
(6.34) 
(6.35) 
(6.36) 
(6.37) 



rid! 


= 2(air 2 + a 2 ri)(cd - b) + 2(rf 


+rf; 


oyrd 


r 2 d' 


= 2{oir 2 +o 2 ri){cd-P)-2{rl 


+ r| 


)a 2 rd 


rid' 


= 2(air 2 + a 2 ri)(cd- (3)-2(r 2 1 


+ r| 


)aird 


r 2 d! 


= 2{air 2 + a 2 r 1 )(cd-b) + 2{r 2 1 


+ rl 


a 2 rd 


T\d 


= 2{oir 2 + o 2 ri)dd + 2o\s 






r 2 c 


= 2{a\r 2 + o~ 2 ri)dd + 2a\s 







and 



(rib + r 2( 3)' 
(n/3 + r 2 b)' 



(nc+ r 2 c)d' - 2(r\ + rj)rai(cd - b) + 2(r\ + r\)ra 2 (cd 
^io 2 +r 2 a 1 ^ 2 ( 
Jir 2 



r i + r 2)d d + {ria 2 + r 2 a{)sd 

ri c + r 2 c)d' + 2(rf + r 2 )aiT(cd- f3) - 2(r\ + r 2 )a 2 T(cd-b) 
ri(J 2 + r 2 ai ( , 2 2 ~ 2 



4cricj2sd (6.38) 



rir 2 



(rf + r\)d 2 d + (ricr 2 + r 2 <7i)a^ - Ao- l0 - 2 sd. (6.39) 



Next we consider the compatibility condition with respect to r. By writing the matrix 
equation (|6.3ip in entrywise form, with the help of Maple, we obtain, with the prime 
denoting the derivative with respect to r, 

-- (rf + ri)(r?TjS + r\r^ + ncjfif + r 2 c(3 + r 2 d 2 d - r lC 2 d + 2sid + r 2 f) (6.40) 



rid' 
r 2 cf 
rid 
r 2 d 
6 



(r l + r 2 )(r l Tb + r 2 r& — ric6 — r 2 cb — rid d + r 2 c d — 2s 2 d — r\f) (6.41) 



(r( + r 2 ){r{rb + r 2 r6 - r 2 cb - r x cb - r 2 d 2 d + r x c l d - 2s x d - r 2 f) (6.42) 



(r'i + r 2 )(r'iTf3 + r 2 r/3 + r 2 c/3 + nc/3 + ncPd - r 2 c 2 d + 2s 2 d + nf) (6.43) 



(r 2 +r 2 )(^-^) 
(rf + rl)(d/3-d6) 



(6.44) 
(6.45) 



and 

ri &' + r 2 /3' 



r 2 6' + n/3' 



(r 2 + r|)(-nc 2 6 + r 2 c 2 f3 - ndd(3 + r 2 ddb + 2si& - 2s 2 6 

- r 2 r/ - r\rf - ncf + r 2 cf) 
(rf + r|)(-r 2 c 2 6 + ric 2 ,S - r 2 ddj3 + r 1 ddb + 2s 2 b- 2s x b 

- rjrf - rfrf - r 2 cf + ncf). 
Direct calculations show that all these equations are satisfied by (|6.4p - (|6.13p . □ 
Lemma 6.7. Theorem 1 6. S\ holds true if r = 0. 

Proof. Equations (|6.4[ )— (|6.7[) follow from [13} Th. 2.4]. The other equations in Theorem l6.2l 
then follow from (^32V(fQ5D and $M$ - §33§ . □ 

With the help of Lemmas I6.6H6.71 one can prove Theorem 16.21 for r ^ in the same 
way as in |14l Sec. 5], where the symmetric case was considered. This is a lengthy and 
tedious calculation that follows exactly the same plan as in |14j . We note that the same 
reasoning also yields the solvability statement in Proposition 12.21 We do not go into the 
details. 
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Alternative approach to Theorem 16.21 

The above reasoning does not give any insight on the origin of the expressions in Theo- 
rem [621 Therefore, in the remaining part of this section, let us deduce these formulas in 
a more direct way. The calculations below are partly heuristic in the sense that we will 
make an ansatz (|6.54p . (|6.62[) . We start with 

Lemma 6.8. The numbers d, d in f)6.1|) satisfy the system of coupled second-order differ- 
ential equations 

d 2 d A . .dd 
— 2 = 4r(noi - r 2 a 2 )- 

- 4(r? + rl){al + al)r 2 d + 8 ^ + ^ d 2 d + S ^ + ^L d, (6.46) 

r±r 2 rir 2 (rf+r|) 



d 2 d , . ,dd 
= -Mnv, - r 2 a 2 )- 



.2 W 2j, (n<T 2 +r 2 ai) 2 ~ 2 ^ g (na 2 + r^) 3 
r\r 2 nr 2 (r 2 +r^) 



- 4(r? + rl){u\ + o\)T l d + 8 v± " z JJ^ iy + 8 ^-^f^ sd. (6.47) 
Moreover, 

3d nr 2 (r 2 + r|) dd 2 2 ^ 2 a\r 2 - o 2 r x 2 

— = ■ r— + (r x + r 2 ) ■ r d + 2{r x s\ - r 2 s 2 )d. (6.48) 

dr cir 2 + cr 2 ri as a±r 2 + cr 2 ri 

Proof. Equation (l6^6|) follows from (102l - (ffn3l ) and (lOHjl - ffOKl) after some lengthy 
algebraic manipulations. Equation (|6,47p follows by symmetry. To obtain (|6.48p . first 
note that (|6.32|) f|6.33j) imply the relations 

r 2 (cd - b) - n(cd -fi) + (r 2 + r\)rd = 0, (6.49) 

dd ~ 
r\r 2 — = (<Tir 2 + o 2 r\){r 2 cd - r 2 b + ricd - n/3) + (r\ + r 2 )(a\r 2 - a 2 ri)rd. (6.50) 

OS 

Now by eliminating / from (|6.4U|) - (|6.41|) we get 

dd ~ ~ 

-77- = (ri+r|)(r 2 6+ri/3)T-(r 2 6-ri/3)(ric+r 2 c)-o!(rfc 2 -r|c 2 )-|-2(riSi-r 2 s 2 )d, (6.51) 

OT 

On account of ()6.49p this becomes 

dd ~ 

q- = ( r i + r I)( r 2 + r iP - r icd - r 2 cd)r + 2(nsi - r 2 s 2 )d. 

Combining this with (|6.50p we obtain (|6.48p . □ 

We seek a solution to the differential equations (|6.46p - (|6.47p in the form 

d = e h g, d = e~ h g, (6.52) 

with h = h(s,r) an odd function of r and g = g(s,r) an even function of r (recall 
Lemma [6. ip . Plugging this in (|6.46p we find, with again the prime denoting the derivative 
with respect to s, 



g" + 2h'g' + {(h'f + h")g = 4r(r 1 a 1 - r 2 a 2 )(g' + h'g) 

2 + r 2 ai 

nr 2 rir 2 (rf+r2 



Af 2 , 2n 2/ 2 . 2^ . Q (n^ 2 + r 2 o-i) 2 3 (ri(J 2 + r 2 CTi) 3 
4(rJ + r 2 )r (<t| + a 2 )g + 8 g a + 8 — ^-sg. (6.53) 
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To obtain further progress we make the ansatz 

d 2 h dh . , 

^2=0, ^ = 2(^-^05)^ (6.54) 

After a little calculation, (|6.53j) then simplifies to 

g = -4(ric72 + r 2 <ri) r g + 8 +8 — — 3— s#. (6.55) 

We can relate (|6.55|) to the Painleve II equation. We have that q = q(s) satisfies q" = 
sq + 2q 3 , if and only if 

g(s) := ciq(c 2 s + c 3 ) (6.56) 



satisfies 



g" = 4c 3 g + 2-| 5 3 + c^ 5 . (6.57) 



Comparing coefficients with (|6.55p . we see that 



(nr 2 ) 1/6 

Cl = (r? + r|)V3 ( 6 - 58 ) 
(ri<r 2 + r 2 (Ti) 
(nr 2 (rf + r^)) 1 ^ 

c 3 = -{rir 2 {rl+rl)f/\ 2 . (6.60) 



Finally, substituting the formulas (|6.52|) . (|6.56j) . 

( nr2 )l/6 



* = *9 = e n ,Y+\ 1/3 « [ 2 ( (r 7 2 2 t r ?W3 S " (™<r? + rl)) 2/3 - 2 ) (6.61) 
in ()6.48p and using again (|6.54p we find after some calculations, 

— = ( r | _ r 4) r 2 + 2 ( riffl _ r2fT2 ) s , (6.62) 

where we are assuming that the choice of the Painleve II solution q in ()6.61|) is independent 
from r. From ()6.61 [) — (|6.62[) and the known result for r = |13t Th. 2.4] we get the 
expression for d in Theorem 12.71 (with q the Hastings-McLeod solution to Painleve II) . 
By symmetry we obtain the expression for d. From (|6.44p . (|6.32p . (|6.34p and a little 
calculation we then find 

JU = -2r- x \rl + rl)TC-\\o) = Ik 1 ^)), 

where the second equality follows from ()4.19p and ()2.27p . Combining this with the known 
result for r = [1 3|. Th. 2.4] we get the expression for c in Theorem 12.71 From (|6.32|) 
and a little calculation we find the expression for b, while (|6.4ip yields the formula for /. 
Finally, the remaining formulas in Theorem 12.71 follow from symmetry considerations, see 
Lemma 16.11 □ 
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